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Abstract 

The profile of a nonlinear stationary thermomagnetic wave in the resistive state of superconductors is 
studied at different transport currents. It is proved that the thermomagnetic wave has an oscillating 
profile at relatively high values of the transport current in the sample. A shock wave with a 
monotonic structure corresponds to comparatively weak transport currents. The wave propagation 
velocity and the wave front width in a superconductor are estimated. 

The problem of thermal suppression of the superconducting state in the presence of transport 
current is topical in the practical use of superconducting magnetic systems (see[l]). External 
fluctuations of various physical natures (thermal, magnetic, mechanical, etc.) can lead to the 
formation of resistive-phase regions in a superconductor. The ohmic heating of these regions 
under a transport current with density jtr can bring about an increase in the temperature 
T above the critical point (T > T^) and, hence, the emergence of a vortex electric field E in 
the sample. It was shown earlier [2] that, together with dispersive and nonlinear effects, the 
dissipative processes induced under ohmic heating in a superconductor lead to the formation 
of stable structures, thermoelectromagnetic E and H waves, depending on the condition at 
the sample surface. In the present work, we studied the qualitative pattern of the dynamics 
of nonlinear thermomagnetic dissipative structures in the resistive state of a superconductor. 
The dynamics of thermomagnetic waves, which move at a constant velocity v along the x 
axis in a superconductor, is described in terms of self-simulated variables ^(x, t) = x — vthy 
the nonlinear heat conduction equation [2-4] 
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and their related equation of the critical state 
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quantity a — 



describes the thermally activated weakening of the Abrikosov vortex 



Here, ly and n are the heat capacity and thermal conductivity coefficients, respectively; 
jc{T) — Jo — a(r — To) is the critical current density in the Bean model — 0^ [5] (where 

To and Tc are the initial and critical temperatures of the superconductor, respectively, and the 

djc_ 

pinning by lattice defects); 2r{E) is the resistive current density; and jtr is the transport 
current density. In the weak-heating approximation (T — Tq) << ^ ~ ^o) for the viscous 
flow region {E > Ef, where Ef is the boundary of the linear section in the current- voltage 
characteristics of hard superconductors [6]), the dependence of the resistive current density 
jr on the electric field E is linear; i.e., jV ~ afE {af is the effective conductivity). In weak 
fields {E < Ef) the dependence jr{E) is essentially nonlinear and is most hkely determined 
by the thermally activated motion of the magnetic flux (flux creep [6]). Here, we will consider 
a perturbation with a sufficiently high amplitude {E > Ef) and use the linear dependence 
jr{E). The appropriate thermal and elecrtodynamic boundary conditions for Eqs. (l)-(4) 
have the form 

T(^ +oo) = To, ^{x -oo) = 0, 

(5) 

E{C ^ +oo) = 0, E{C ^ -oo) = E„ 

where E^ is the external electric field. 

By solving the set of Eqs. (l)-(4) with boundary conditions (5), we obtain the following 
equation for the ditribution of the nonlinear E wave: 
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Here, z = —,/3 = ^^,Ej = — , r = — - — = — 5-, L = — — - is the depth of the 
L L af c^u aL^ 47rjo 

magnetic field penetration into the bulk of the superconductor, — is the thermal 

K 

time of the problem, and Hg is the external magnetic field. 
By analyzing the roots of the characteritic equation [7] 

X{X + P){X + Pt)^-P't^, (7) 

we can easily verify that the problem considered has solution that oscillate about the origin 
of the coordinates {Eq = 0) of the phase diagram with decreasing and increasing amplitudes. 

The oscillations appear at Ej > — rE^, attenuate at 3(1 + t)Ek > E,; > — tT'k, 

J 2(1 + r) i-^y I ^ J 2(1 + r) ' 

and increase at f3{l + t)Ei^ < Ej. This means that, at a sufficiently high transport current 
itr-i the roots of the characteristic equation (7) are complex, the integral curves are helical. 
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and the origin £'0 = is a singular point (focus) . If the coefficients (5 and r are assumed to 
be universally positive, the representative point approaches the origin of coordinates, which 
is a stable equilibrium position. In the opposite case, when j3{l -\- t)Ek < Ej, the focus 
becomes unstable and the representative point goes to infinity. Note that the final value of 
the transport current jtr leads to a displacement of singular points and, correspondingly, of 
the boundary conditions. 

dE 

At sufficiently weak transport currents, when the condition — — >> P'^tEj is satisfied 
(high-amplitude waves), Eq. (6) is an integrable equation and the corresponding model can 
be considered a damped linear oscillator owing to the friction force, that is, 

F, = -/3(l + r)^, (8) 

where z is an analogue of time and E is the coordinate of the "material" point. The potential 
well equation has the form 

dE 

Analysis of the phase plane ("^) E) described by Eq. (6) shows that it has two equilibrium 
points: Eq is a stable node and Ei = 20^tEi^ is a saddle. The two equilibrium states {Eq, Ei) 
are separated by the scparatrix AB in the phase plane (Fig.l). The material point is located 
at the point £"1 at 2; —> 00 and goes over into the point Eq ai z ^ —00. The change-over 
from one equilibrium state to other occurs only monotonically. 

Within the approximation r << 1, the solution of Eq. (6) is represented as 



1 - th{^{z - Zq) 



(10) 



The condition r = — — << 1 means that the magnetic flux is redistributed much faster 
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than the heat is transferred. Here, Dt — — and — are the thermal and magnetic 

diffusion coefficients, respectivily [6]. Hence, the spatial scale Le for the magnetic flux 
penetration is substantially greater than the corresponding thermal scale Lt- Therefore, the 

spatial derivatives -^--^ will contain a small parameter 

validity of this approximation by the direct differentiation of Eq. (6), that is, 

d'E 



<< 1. We can easily check the 




(11) 



The maximum error of this approximation is of the order ^ — — . For example, the error 

is equal to 25 small in the hmiting case r — > 0. Relationship (10) describes the profile of 
a shock thermomagnetic wave propogating into the bulk of the superconductor. The wave 
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structure is shown in Fig. 2. At the boundary condition {E{z — > — oo) = Eg), we can easily 
determine the velocity ve oi the wave with an amplitude Eg in the form 

L 
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(12) 



and the wave front width 5z 



Sz = 16 
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The numerical estimates yield the values ve — 1 10^ and 5z — 10~^ 10~^ for r — 1. 

In conclusion, we note that the results obtained make it possible to describe the nonlinear 
stage of evolution of the thermomagnetic instability in the resistive state of superconductors. 
In the initial stage, wc observe an exponential increase in the perturbations of T, E, and 
H with the increment A determined from the linear theory (see [8]). The linear stage of 

the instabihty evolution lasts for the period tj — —. In the latter stage of the instability 

— * — * — * — * 

evolution, the stationary E or H wave (depending on whether the E or H quantity is fixed 
on the sample surface during the wave motion) propagates into the bulk of the sample. The 
velocity of the wave is determined by formula (12), and the time of its motion in the sample 
with thickness 21 is given by 

Ve 

Hence, it follows that the transition from the superconducting to the normal state can occur 
through the propagation of a stationary thermomagnetic wave whose structure essentially 
depends on dispersive and dissipative effects. 
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